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Abstract 



The Collins and Sivers asymmetries for charged hadrons produced in deeply inelastic 
scattering on transversely polarised protons have been extracted from the data 
OJ ■ collected in 2007 with the CERN SPS muon beam tuned at 160 GeV/c. At large 

o,. values of the Bjorken x variable non-zero Collins asymmetries are observed both for 

O I positive and negative hadrons while the Sivers asymmetry for positive hadrons is 

slightly positive over almost all the measured x range. These results nicely support 
the present theoretical interpretation of these asymmetries, in terms of leading- twist 
quark distribution and fragmentation functions. 
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After first indications of transverse spin effects in liadron physics in the 1970s [H |2] 
their importance was unambiguously estabhshed by the remarkably large single spin asym- 
metries (SSAs) found in pp collisions at Fermilab both for neutral and charged pions [3]. 
Following the discovery by the EMC at CERN in 1988 that the quark spins contribute 
only little to the proton spin |1] , the interest in the nucleon spin structure was revived and 
a more complete description including quark transverse spin and transverse momentum 
has been worked out. 

The quark structure of the nucleon in the collinear approximation or after integra- 
tion over the intrinsic quark transverse momentum kx is fully specified at the twist-two 
level by three parton distribution functions (PDFs) for each quark flavour j[5j : the momen- 
tum distributions q{x), the hehcity distributions Aq{x) and the transverse spin distribu- 
tions Axqix), where x is the Bjorken variable. The latter distribution — often referred to 
as transversity — is chiral-odd and thus not directly observable in deep inelastic scattering 
(DIS). In 1993 it was suggested \6\ that transversity could be measured in semi- inclusive 
lepton-nucleon scattering (SIDIS) due to a mechanism involving another chiral-odd func- 
tion in the hadronisation, known today as the Collins fragmentation function (FF). The 
mechanism leads to a left-right asymmetry in the distribution of the hadrons produced 
in the fragmentation of transversely polarised quarks. Thus a transverse spin dependence 
in the azimuthal distributions of the final state hadrons can be generated both in trans- 
versely polarised pp scattering and in SIDIS off transversely polarised nucleons. In the 
latter case the measurable Collins asymmetry, Acoii, is proportional to the convolution of 
the transversity PDF and the Collins FF. 

Admitting a finite kx, in total eight PDFs are needed for a full description at leading 
twist and leading order in a^ [3 El E] ■ All these functions lead to azimuthal asymmetries 
in the distribution of hadrons produced in SIDIS processes and can be disentangled mea- 
suring the different angular modulations. Amongst the transverse momentum dependent 
PDFs, the T-odd Sivers function [10] is of particular interest. This function arises from a 
correlation between the transverse momentum of an unpolarised quark in a transversely 
polarised nucleon and the nucleon spin. It can be different from zero because of final state 
interactions mediated by soft gluon exchange between the interacting quark and the tar- 
get remnants [TT]. It is responsible for the Sivers asymmetry, Asiv, which is proportional 
to the convolution of the Sivers function and the unpolarised FF. The Sivers mechanism 
might also be the reason for the large asymmetries observed in pp collisions. 

Transverse spin effects in SIDIS are investigated, at different beam energies, by the 
HERMES experiment at DESY and the COMPASS experiment at CERN. An experiment 
to measure transversity using a transversely polarised ^He target has recently been per- 
formed at JLab [12j. Transverse spin effects are also an important part of the scientific 
programme of the RHIC spin experiments at BNL. 

Up to now, sizable Collins asymmetries for the proton were observed recently by 
HERMES using a proton target [13]. This implies non- vanishing Collins fragmentation 
and transversity functions. Direct measurements at the KEK e~^e~ collider by the BELLE 
experiment established that this Collins FF is sizable [lll[T5]. COMPASS measured van- 
ishing asymmetries by scattering high energy muons off a deuteron target [161 [13, HI]. All 
these data were well described by a global fit [191 [20] which allowed for a first extraction 
of the u and d-quark transversity PDFs. 

The Sivers asymmetry for the proton was measured by HERMES [T31 UT] to be 
different from zero for positive hadrons, while it was found to be compatible with zero 
for deuteron by COMPASS [IH [13 [E] . These HERMES and COMPASS data could also 
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be well described by theoretical calculations and fits, and allowed for extractions of the 
Sivers function |22], which turned out to be different from zero and opposite in sign for u 
and (i-quarks. 

In this Letter, we present the COMPASS results on the Collins and Sivers asym- 
metries for charged hadrons produced in SIDIS of high energy muons on transversely 
polarised protons. The data were collected in 2007 using NH3 as target material and a 
160 GeV/c beam with a momentum spread Ap/p = ±5%. The beam was naturally po- 
larised by the vr-decay mechanism, with a longitudinal polarisation of about -80%. This 
measurement followed the measurements performed in 2002, 2003 and 2004 at the same 
energy with the transversely polarised ^LiD target. 

The COMPASS spectrometer [23] is in operation on the M2 beam line of CERN since 
2002. Two magnetic stages are used to ensure large angular and momentum acceptance. 
A variety of tracking detectors is used to cope with the different requirements of position 
accuracy and rate capability at different angles. Particle identification is provided by a 
large acceptance RICH detector, calorimeters, and muon filters. Major upgrades in 2005 
mainly concerned the polarised target, the tracking system, the RICH detector, and the 
electromagnetic calorimeters. The new target solenoid magnet provides a field of 2.5 T 
and has a polar angle acceptance of 180 mrad as seen from the upstream end of the target. 
In the earlier measurements with the ^LiD target the polar angle acceptance was 70 mrad. 
The target material is cooled in a ^He-^'He dilution refrigerator, and the protons in the 
H atoms are polarised to 0.80-0.90 by dynamical nuclear polarisation. About 48 hours 
are necessary to reach 95% of the maximal polarisation. A pair of saddle-shaped coils 
can provide a 0.6 T vertical field which is used to rotate the target nucleon spin and to 
hold the polarisation vertical for the transversity measurements. In the frozen spin mode, 
and with the holding field at its operational value, the relaxation time of the polarisation 
exceeds 3000 hours. 

The target consisted of three cylindrical cells with 4 cm diameter, one central cell 
of 60 cm length and two outer ones of 30 cm length, all separated by 5 cm. Neighbouring 
cells were polarised in opposite directions, so that data with both spin directions were 
recorded at the same time. In order to minimise the effects due to different spectrometer 
acceptance for different target cells, in each period of data taking a polarisation reversal 
was performed after 4-5 days by changing the microwave frequencies in the three cells. 

The geometry of the polarised target and the data taking procedure were chosen 
such as to optimise the extraction of spin asymmetries. The principle of the measurement 
can be understood by considering the "ratio product" [T7] 

p inner outer /-, \ 

where Nj^^^^ and Neuter are the number of hadrons produced in the first sub-period on 
oppositely polarised cells, and A^^tiner ^^^ A'oMier ^^^ the corresponding numbers in the 
second sub-period, i.e. after polarisation reversal. The ratio product is constructed such 
that beam flux, spin-averaged cross-section, and the number of scattering centres cancel. 
As long as the ratios between the spectrometer acceptances of each cell are the same in 
the two sub-periods and the number of produced hadrons follows the generic azimuthal 
modulation N"^^ ~ 1 ± e sin $, one simply gets i? = 1 + 4e sin $, and the extraction of the 
amplitude e of the azimuthal modulation is straightforward. 

In 2007 data were taken at a mean beam intensity of about 5 x lO'^ /^"""/s (typically 



2.4 X 10^ /i"^/spill, for a spill length of 4.8 s every 16.8 s). Using up 4 x 10^'^ niuons, about 
12 X 10^ events were collected in six separate periods, corresponding to 440 TB of data. 

In the data analysis, events were selected if they had at least one "primary vertex" , 
defined as the intersection point of a beam track, the scattered muon track, and other 
possible outgoing tracks. The momenta of both incoming and outgoing charged particles 
were measured. The primary vertex was required to be inside a target cell. In order to 
guarantee the same muon flux along the target material, the extrapolated beam track 
had to traverse all the three target cells. For incoming and scattered muon tracks, as well 
as for the other reconstructed tracks, y^ cuts were applied to assure the quality of track 
reconstruction. Tracks from the primary vertex which traversed more than 30 radiation 
lengths were identified as scattered muons. The event was rejected if more than one of 
such tracks were found. 

In order to be in the DIS regime, only events with a photon virtuality Q^ > 
1 (GeV/c)^, a fractional energy of the virtual photon 0.1 < y < 0.9, and a mass of 
the hadronic final state W > 5 GeV/c^ were considered. The variable x covers the range 
from 0.004 to 0.7. 

All particles emerging from the primary vertex were assumed to be hadrons if they 
traversed less than 10 radiation lengths of material. For tracks with an associated cluster 
in one of the hadronic calorimeters, a minimal amount of deposited energy was required 
to further reduce the electron and muon contamination. Finally, tracks reconstructed only 
in the fringe field of the first analysing magnet of the spectrometer were rejected. This 
roughly corresponds to a cut at 1.5 GeV/c in the hadron momenta. In order to recon- 
struct the hadron azimuthal angle with good precision, the hadron transverse momentum 
with respect to the virtual photon direction, p^, was required to be above 0.1 GeV/c. A 
minimum value of 0.2 for z, the relative energy of the hadron with respect to the virtual 
photon energy, was chosen to avoid hadrons from the target fragmentation region. 

As explained in detail in Ref. [TTj, the Collins effect shows up as a modulation 
[1 + ec sm{(j)h + 05 — vr)] in the number of events, where (ph and 4>s are the azimuthal 
angles of the hadron and of the target nucleon spin vector in a reference system in which 
the z-axis is the virtual photon direction and the x-z plane is the lepton plane according 
to Ref. [21]. . The amplitude of the modulation is ec = DmnIPtA-coIU where -Dattv = 
(1 — y)/(l — y + y^ 1'^) is the transverse spin transfer coefficient from target quark to 
struck quark, / the dilution factor of the NH3 material, and Pt is the proton polarisation. 
Similarly, the Sivers effect results in a modulation [l+es sin(0;i— 05)], where e^ = fFxAsiv 

The transverse spin asymmetries were obtained by comparing the azimuthal distri- 
butions of the detected hadrons as measured in the first sub-period of data taking with 
the corresponding distributions of the second half measured with opposite target polari- 
sation. Since the two sets of data were taken typically one week apart, the stability of the 
apparatus is a central point in the measurement. As a first step in the data selection, the 
hit distributions of all trackers were scrutinised, as well as the number of reconstructed 
events, the number of vertices per events, and the number of tracks per event. In a second 
step, the stability of the average tt+tt" invariant mass in the i^° region as well as the 
distribution of twelve kinematic quantities (x, y, W, z, ...) were investigated dividing the 
data in small time-ordered sub-samples. Each distribution of each sub-sample was com- 
pared with the corresponding ones of each other sub-sample within the same data taking 
period, and sub-samples were rejected when deviating more than 3.5 agtat from the mean 
values. 

As a final selection criterion, the data were tested for a possible dependence on 



either sin(0/i + (ps) or sin^cph — (ps) of the acceptance ratio between two consecutive sub- 
periods with opposite target polarisation. Combining the number of events reconstructed 
in the different target cells in two consecutive data taking sub-periods, one can construct 
two different estimators on the stability of the acceptance. The first estimator measures 
the mean modulation in the relevant azimuthal angle of the acceptance ratio between two 
sub-periods. The second one probes possible large differences in the acceptance ratios for 
the different target cells which could affect the physics asymmetry. These two pieces of 
information have been used to construct a x^ ^^^ fhe final selection of the data taking 
periods was done on the basis of its value. 

As a result of the quality control, all data collected in the six periods were used for 
the extraction of the Collins asymmetry, while only four periods were used for the Sivers 
asymmetry. This can be understood because the Sivers asymmetry is very sensitive to 
instabilities of the spectrometer since it is the amplitude of a modulation of the azimuthal 
angle of the hadron transverse momentum with respect to the target spin vector. On the 
contrary, the Collins asymmetry is an asymmetry in the azimuthal angle between the 
hadron transverse momentum and a direction which depends on the target spin direction 
and the lepton scattering plane, which is different for each event. The final sample contains 
23.1 X 10^ SIDIS events for the Collins asymmetry and 15.6 x 10^ for the Sivers asymmetry. 

The asymmetries were evaluated for positive and negative hadrons in bins of the 
three kinematic variables x, z and j9y. The binning is the same as used for the previous 
analyses of deuteron data and consists of 9 bins in x, 8 bins in z and 9 bins in p^, 
integrating over the other two variables. For each period, the physics asymmetries were 
obtained by dividing the raw asymmetries by the target polarisation, the dilution factor, 
and, in the case of Collins analysis, by the -Datat factor. The target polarisation was 
measured individually for each cell and each period. The dilution factor of the ammonia 
target was evaluated for each bin. It is 0.15 in average, and increases with x from 0.14 to 
0.17. 

The estimator used for the evaluation of the raw asymmetries is based on an ex- 
tended unbinned maximum likelihood method [25j. The likelihood function is built as the 
product of the probability densities p corresponding to each hadron i from each target 
cell. The likelihood for hadrons from a given target cell in one period is written as 

(N+ \ ivT / Af- 

i=0 J \ 4=0 

The + and — signs refer to the orientation of the target polarisation in the two sub- 
periods and A^"*^ is the corresponding total number of hadrons. The quantities J^ are the 
integrals of the probability densities over (ps and (ph- The probability densities p"^ are 
the product of two parts, one corresponding to the acceptance description and the other 
to the SIDIS cross section of longitudinally polarised leptons on transversely polarised 
nucleons. Various parametrisations of the acceptance part were tested, resulting in a 
negligible dependence of the extracted asymmetries on the acceptance description. The 
cross section was parametrised taking into account both the unpolarised and polarised 
parts. The polarised part consists of all the expected eight modulations, namely sm{(ph + 
(ps-rc), sm{(ph - (ps), cos{(ph - (ps), sin(20/i - (ps), cos{2(ph - (ps), sin(05), cos(05), and 
sin{3(ph — (ps)i aiid all their amplitudes were extracted at the same time. The Collins and 
Sivers asymmetries are proportional to the amplitudes of the first two terms. Systematic 
studies for the other six amplitudes are still ongoing, and those results will be the subject 



"^^ 

^ 



0.1- 



i\ 



^ 



*» 



• positive hadrons 
O negative hadrons 



-^#^^ 



^ 



^^^ 



-0.1 



[=:=: 



M 



ri 



_L 



10" 



10" 



0.5 



0.5 



P 



1.5 
'' (GeV/c) 



Figure 1: Collins asymmetry as a function of x, z, and p^, for positive (closed points) and 
negative (open points) hadrons. The bars show the statistical errors. The point to point 
systematic uncertainties have been estimated to be 0.5 o-gtat for positive and 0.6 astat for 
negative hadrons and are given by the bands. 



of a future publication. 

The final Collins and Sivers asymmetries extracted with the likelihood method were 
compared with the asymmetries extracted using four other estimators, including those 
used in the previous publications which were based on the "ratio product" R of Eq. [H 
finding an excellent agreement between all results. The correlation coefficient between the 
Collins and the Sivers asymmetries turned out to be small, less than 0.2 in absolute value 
over the whole x range. 

Extensive studies were performed in order to assess the systematic uncertainty of 
the measured asymmetries. All the studies were done separately for positive and negative 
hadrons and for the Collins and the Sivers asymmetries. 

The largest systematic error is due to residual acceptance variations within pairs of 
data taking sub-periods. To quantify these effects, two different types of false asymmetries 
were calculated, using the external cells and the internal cell divided in two parts, and 
assuming wrong sign polarisation for one of the two. Moreover, the physical asymmetries 
were also extracted using only the first and only the second half of the target. The dif- 
ference between these two physical asymmetries, the false asymmetries, and the degree of 
compatibility of the results from different periods were all used to quantify the systematic 
uncertainty. 

In the case of the Collins asymmetry, the systematic uncertainty is estimated to be 
0.5 astat for positive and 0.6 astat for negative hadrons. In the case of the Sivers asymmetry, 
the systematic error is 0.8 astat for positive and 0.4 astat for negative hadrons. A further 
systematic uncertainty of ±0.01 is present in the absolute scale of the Sivers asymmetry 
for positive hadrons. It reflects a 0.02 difference in the mean value of the asymmetries 
extracted in the first two and in the second two periods of data taking used for this 
analysis. In spite of throughout studies, the origin of this difference, which affects only 
the Sivers asymmetry for positive hadrons, could not be identified and had therefore to 
be included in the systematic uncertainty. The results of this measurement of the Collins 
and Sivers asymmetries are shown in Fig. [T] and [2] as a function of x, z, and p^, for 
positive and negative hadrons. Figure [3] displays the mean values of kinematic variables 
for positive hadrons in the x, z, and p^ bins. The corresponding quantities for negative 
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Figure 2: Sivers asymmetry as a function of x, z, and p^"? for positive (closed points) and 
negative (open points) hadrons. The bars show the statistical errors. The point to point 
systematic uncertainties have been estimated to be 0.8 o-gtat for positive and 0.4 astat for 
negative hadrons and are given by the bands. For positive hadrons only, an absolute scale 
uncertainty of ±0.01 has also to be taken into account. 



5 25 

> 

D on 

O 20 


- T <p';>(GeV/c) 
O (z) 

- • <e^>(GeV/c)- • 


- ▼ <p'')(GeV/c) 
(xfx 10 
• (e'XGeV/cf 




- ▼ (x) X 10 
(z) 

- • (e^XGeV/c)^ 


~ 


C 15 
^ 10 


" 000000080 


y T ▼ 

-(8553 


T 



"00000 


9 - 


5 


• 
• 

• f • •* 

r 1 


-•••• • • • 




• 
1 1 


••••• • • • 

1 1 


• 



10" 



10-' 



0.5 



0.5 



0.8 O 



0.6 -=a,^ 
0.4 ^ 



0.2 



1.5 
>'' (GeV/c) 



Figure 3: Mean values of some kinematic variables in the final data sample. From left to 
right: mean values of p^, z and Q^ as functions of x; mean values of p^, x and Q^ as 
functions of z; mean values of x, z and Q^ as functions of p^. 



hadrons are very similao. 

As it is clear from Fig. [H the Collins asymmetry has a strong x dependence. It 
is compatible with zero at small x within the small statistical errors and increases in 
absolute value up to about 0.1 for x > 0.1. There, the values agree both in magnitude 
and in sign with the previous measurements of HERMES [13], which were performed at 
the considerably lower electron beam energy of 27.5 GeV. Also, the present results agree 
with the predictions of the global analysis of ref. [HI |20] and thus strongly support the 
underlying interpretation of the Collins asymmetry in terms of a convolution of the twist- 
two transversity PDF and the FF of a transversely polarised quark. An important issue 
is the Q^ dependence of these functions. Our results at large x are compatible with the 
HERMES data in spite of the higher Q^ values which exceed those of HERMES by a 
factor 2 to 3 with increasing x. This indicates that the possible Q^ dependence should not 



All numerical values have been put to HEPDATA. 
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Figure 4: Collins (upper row) and Sivers (lower row) asymmetry as a function of W, for 
positive (left) and negative (right) hadrons. The closed and open points give the values 
for the "large x" and the "small x" samples respectively. The errors are statistical only. 



be dramatic in the present energy ranges. 
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Figure 5: Mean values of Q^ (left) and x (right) as functions of W. The closed and open 
points give the values for the "large x" and the "small x" samples respectively. 



The results for the Sivers asymmetry for negative hadrons exhibit values compatible 
with zero within the statistical accuracy of the measurement. For positive hadrons, the 
data indicate small positive values, up to about 3% in the valence region. These values 
are somewhat smaller than but still compatible with the ones measured by HERMES at 
smaller Q^. Given the importance of the Sivers function in the present description of the 
transverse momentum structure of the nucleon, we looked at a possible kinematic depen- 
dence of our measurements. In particular, we evaluated the asymmetries as a function of 



W, separately for the "large-x" (x > 0.032) and "small-x" (x < 0.032) samples. The re- 
sults are shown in Fig. HI The mean values of Q^ and x in all W bins are given in Fig. [51 As 
it is apparent from Fig. HJ no conclusion can be drawn about a possible W dependence of 
the Collins asymmetry. On the other hand, the signal of the Sivers asymmetry for positive 
hadrons seems to be concentrated at small W, in the region where HERMES measures, 
and goes to zero at large W, which for large x means large Q^. Thus our data give an 
indication for a possible W dependence of the Sivers asymmetry for positive hadrons. 
Definite conclusions will be possible only when new more precise data at high energy will 
become available. 

In summary, for the first time the Collins and Sivers asymmetries for positive and 
negative hadron production in DIS off the proton have been measured at high energy. 
Our data extend the kinematic range to large Q^ and large W values. The x range has 
been extended to considerably smaller values which are needed to evaluate the PDF first 
moments. For the Sivers asymmetry, a signal is seen for positive hadrons, which persists 
to rather small x values. The data give an indication for a possible W dependence of this 
asymmetry, but the present statistical and systematic uncertainties do not allow definite 
conclusions. The measured Collins asymmetry is sizable for both positive and negative 
hadrons also at high energies and Q^. Thus Collins asymmetries measured in SIDIS are 
an appropriate tool to investigate the transversity PDF. 
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